**Volume Title** 

ASP Conference Series, Vol. **Volume Number** 
**Author** 

© **Copyright Year** Astronomical Society of the Pacific 



Clumpology of Starbursts in the WiggleZ Dark Energy Survey 

Emily Wisnioski, Karl Glazebrook, Chris Blake 

Centre for Astrophysics and Supercomputing, Swinburne University of 
Technology, P.O. Box 218, Hawthorn, VIC 3122, Australia 
( ewisnios @ astro, swin. edu. au ) 

& The WiggleZ Team 

Abstract. We have observed Ha emission from a sample of 13 star-forming galaxies 
from the WiggleZ Dark Energy Survey in the redshift range z ~ 1.2 — 1.5 and stellar 
mass range 9.8 < log(M„/M )< 11.6 with the integral field spectrograph OSIRIS on 
Keck, taking advantage of Laser Guide Star Adaptive Optics. We detect multiple emis- 
sion, 1-2 kpc size sub-components, or 'clumps' within the Ha spatial emission in 4 
galaxies, extended emission in 4 galaxies and compact spheroidal emission in 5 galax- 
ies. Considering our results in the context of unstable disk formation and a merger 
sequence, we find evidence of ordered orbital motion in the majority of galaxies as 
would be found in unstable gaseous disks. 



1. Morphologies at High Redshift 

Morphological studies made possible by deep imaging surveys have provided the first 
clues to the processes of high-redshift galaxy formation. It has become clear from these 
results that high-redshift galaxies exhibit different morphologi es to local galaxies, with 



higher fra ctions of b oth compact and irregular morphologies (IGlazebrook et al.lll995 
Abra ham et al.lll996h . Initial interpretations of the complex structures seen in the Hub 
ble Ultra Deep Field favor higher rates of major merg ers, a natural consequence of hi 



erarch ical merging models in a ACDM Universe (e.g. JConselice et ai1l2003l : lEotz et al. 



2006), however, additional obse rvations of irregular galaxi es reveal properties reminis- 
cent of local disk galaxies dElmegreen & Elmegreenll2005l) . Yet, unlike local disks, the 
galaxies at high redshift contain large dense re gions or 'clumps' of star formation, com- 
parable to HII regions but of much larger size ( Jones et ai1l2010h . Models have been de- 
veloped to explain the clumpy str uctures in disks and their evolution into galaxies with 
morphologies seen at z = (e.g., lNoguchilll999l : llmmeli et al. I2004al : Bournaud et al. 



120071 : lElmegreen et al.ll2008h . 

Morphological studies alone do not provide enough information to distinguish 
between these kinematically distinct theories of galaxy formation. Spatially resolved 
spectroscopy is therefore essential for disentangling kinematics and star formation his- 
tories of these complex structures at high redshift. 
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2. Kiloparsec scale Observations of Star-Forming Galaxies at z > 1 



The past five years have seen a great number of integral field spectro scopy (IFS) studies 
of nebular emission line Ha in z > 1.5 star-forming ga l axies (e.g., Genzel et al. 2006; 
LawetalJl2009l : IWrieht et alJl2009t lEpinat et all 120091 : iNesvadba et al.ll2008h . These 
redshifts are favored for the large number of Lyman break galaxies (LBGs) and sub- 
millimeter galaxies (SMGs) being discovered by methods tailored for high redshift and 
for the accessibility of Ha in the near infrared. 

A main focus of IFS studies has been the classification of galaxies kinematically to 
define clear observ ational distinctions be t ween disks and in teracting systems. Mergers 
observed with IFS (INesvadba et al.ll2008c IWright et al.ll2009f) have been found to show 
clear kinematic steps across one resolution element indicating a velocity separation be- 
tween two kinematically separate systems. Some galax ies also exhibit rotational signa- 
tures across IFS results and are well fit by disk models (Ge nzel et al.ll2006l : lEpinat etal 



20091 : Lemoine-Busserolle & Lamareiilell2010l) . Yet rigidly classifying galaxies in these 
two distinct groups is an over-simplification. Many of the samples of > 5 galaxies 
present a variety of kinematic types including dispersion dominated systems, stable 
disks, unstable disks, minor mergers, merger, and merger remnants with definitions of 
each class varying from one sample to another. 



3. Discovering extreme objects at z > 1 using the WiggleZ Dark Energy Survey 

We present integral field spectroscopy of Ha emission in 13 star-forming galaxie s se- 
lected from the WiggleZ Dark Energy Survey at z ~ 1.3 (Drinkwa ter et al.ll2010h . By 
moving to slightly lower redshifts than previous studies of LBGs and SMGs we are able 
to study the internal kinematics of galaxies of comparable intrinsic luminosities (Fig- 
ure [B but achieving higher spatial resolution and comparable signal to noise in a much 
shorter exposure time. WiggleZ galaxies are UV selected with color cuts to select the 
bluest galaxies across the whole redshift range of the survey. While the WiggleZ red- 
shift distribution peaks at z ~ 0.6, there exists a long thin tail of high-redshift galaxies 
that reaches as high as z = 1.5. These galaxies are identified from strong [Oil] emission 
and have star formation rates of 50-500 M Q yr _1 . There are about 100 such objects in 
the survey. The WiggleZ survey has made the current study of intensely star-forming 
galaxies at 1 < z < 2 possible by finding rare starbursts that would otherwise be missed. 
These objects may be the link needed to bridge the gap between observations of local 
LIRGs/ULIRGs and high-redshift LBGs and SMGs at one of the most important epochs 
of mass assembly. 



4. Results 

We observed 13 star-forming galaxies using OSIRIS (OH Suppressing InfraRed Imag- 
ing Spectrograph) with the Keck II laser guide star adaptive optics (LGSAO) system. 
By using an IFU we obtain an Ha spectrum for every spatial pixel (spaxel) of the galaxy 
image. A Gaussian profile is fitted to the Ha emission in each spaxel of the cube to map 
out the Ha flux distribution of the galaxy. A velocity map is constructed for each galaxy 
using the velocity shifts of the emission in each spaxel relative to the galaxy redshift. 
Velocity dispersion is measured in each spaxel from the r.m.s. of the fitted Gaussian 
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Figure 1. Ha luminosity as a function of redshift for the largest high-redshift 
IFS surveys with our work shown by the red diamonds. The gray points show the 
[Oil] luminosity of o ther WiggleZ galaxi e s considered for t h is sample. The remain- 
ing symbols are from l Epinat et al.l (2009); L aw et al.1 (120091) : iForster Schreiber et ail 
(2009); Wright etal] (120091) . 



line profile corrected for instrumental broadening (R=3600). The kinematic maps are 
shown for three galaxies in Figure |2] 

We classify the 13 galaxies into three classes based solely on their morphologies; 
single emission: showing a single resolved source of Ha emission, extended emission: 
showing multiple unresolved or extended components, and multiple emission: having 
multiple resolved regions of Ha emission. These classifications are used to discuss 
the properties of galaxies with like morphologies in the context of different formation 
scenarios. 

The Clumpy Disk Hypothesis: We investigat e the scenario that the c lasses repre- 
sent different stages of the clumpy disk models of lElmegreen et all (120081) simulations 
(EBE08). In these simulations gravitational instabilities result in multiple star-forming 
complexes forming under Jeans collapse within a young gas rich disk. As turbulent 
speeds decrease relative to the rotational speed of the disk the clumps migrate to- 
wards the center on a timescal e of ~ 1 Gyr. Finally, the clumps coalesce and form 
a slo wly rotating p seudo -bulge (Noguchi 1999 ^ Bournaud et ai1l2007l : TElmegreen et al. 
20081: iDekel etal1l2009l : iKrumholz & Dekelll2010h . 



The morphology of multiple emission systems is well matched to the early stages 
of clumpy disk formation (t - 175 - 450 Myr) of the EBE08 simulations. The velocity 
of the clumps follow s the global velocity structure of the galaxy, consistent with models 
(llmmeli et al.ll2004bh . This is seen in WK0909_06R (Figure [2]), where the major axis of 
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Figure 2. OSIRIS Ha kinematics of 3 out of 13 galaxies in our sample, 
WK0909_06R is an example of the multiple emission galaxies, WK0912_01R is an 
example of the extended emission galaxies, and WK1002_42S is an example of the 
single emission galaxies. Left: Ha flux map, middle velocity map from the systemic 
velocity, right velocity dispersion map. The black crosses mark the peak(s) in Ha 
flux. Velocity and velocity dispersion are measured in km/s. 



rotation cuts across multiple clumps revealing no corresponding deviation in the rota- 
tion profile. The morphologies of the extended emission galaxies are either consistent 
with the system being observed at the clump coalescence phase (t = 600 - 800 Myr) 
or the result of projection effects, with some galaxies well fit by disk models. Compact 
single emission galaxies have morphologies well matched to the later stage of simulated 
formation when the clumps have collapsed to form a bulge (? = 1000 Myr). The high 
dispersions and masses of these systems are consistent with EBE08 bulges, predicted 
to become slow rotators with high turbulence of <t\ os ~ 130 km s -1 in the inner regions 
with M* > 10 10 M . 
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Star-forming regions within a disk should have similar or consistent metallici- 
ties as they are evolving within the same conditions (llmmeli etal.l f2004a) In a ma- 
jor or minor merger no similar correlation is expected in the metallicities of each 
galaxy. By integrating the light of each clump the clump metallicities are calculated 
from the Ha/[NII] ratio. Metallicities of clumps in multiple emission galaxies, such as 
WK0909_06R (Figure [2]) are equal within the errors, as predicted for the early stages of 
disk fragmentation. 

The Merger Hypothesis: However, with separations of 3-10 kpc the Ha com- 
plexes may be sites of multiple regions of star formation merging from different par- 
ent galaxies. Without deep imaging of diffuse continuum light surrounding the Ha 
emission to reveal tidal features characteristic of mergers, we need to utilize kinematic 
quantities to find support for this hypothesis. We analyze here a scenario in which our 
sample form a merger sequence with pre, ongoing, and post mergers. In this scenario 
the multiple emission galaxies are the initial stages of mergers when star formation in 
separate galaxies can be spatially differentiated. The extended emission galaxies are 
the next stage in the sequence representing ongoing merging activity (or pre-mergers 
with component confusion due to projection effects). The final stage in the proposed 
merger sequence is the single emission galaxies which represent merger remnants with 
high dispersions and star formation rates. 

Although some of the galaxies in our sample show indications of ordered rotation 
indicative of disks from their velocity maps, none of the systems show a clear peak 
in velocity dispersion corresponding to the turnover of the velocity distribution as ex- 
pected for disk galaxies. Although this may be explained by sensitivity effects, the 
lack of a peak in velo city dispersion, a feat ure observed in IFS studies in both local 
( de Zeeuw et al.l l2002h and high redshift dGenzel et al.1 120061 : iForster S chreibe r et al. 



2009) disk galaxies, may indicate that we are seeing separate kinematic components 



merging. In this case, the observed peaks in dispersion may correspond to the centers 
of individual galaxies merging or be a result of shock fronts between merging systems. 
The majority of IFS results that have been classified as mergers show velocity maps 
where the velocity does not vary smoothly a cross the source but r ather show a quick 
turnover within a resolution element or less dNesvadba et al .1120081 : IWright et al.ll2009l : 
Law et al1l2009l) . This step function is not present in the kinematics in our sample. 



Findings: Whilst it is natural to attempt to identify kinematic observations as 
disks or mergers, this simplified picture of galaxies becomes increasingly fallible at 
high-redshift. We are hesitant to split our data as has been done in previous studies 
(disk, merger, dispersion dominated, etc.) on an individual basis as we find that galax- 
ies do not populate homogeneous classes but rather form a heterogeneous sample. In a 
sample of 13 galaxies it is likely a variety of mechanisms are at work, even if disk insta- 
bility is a dominant mechanism in star formation histories at high redshift, as our data 
supports, the disks are likely to undergo mergers in their lifetime further complicating 
kinematic interpretations. Future deep optical and infrared imaging will help to clar- 
ify the extended morphologies of these systems and assist in discriminating between 
turbulent disks and ongoing major/minor mergers. 
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